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Abstract 
> 

O ' The deterministic single-photon emission by means of STIRAP through 

'"T' , the atoms with degenerate levels is studied. The expression for the polar- 

ization matrix of the emitted photon is obtained and its dependence on 
the polarization of the driving laser field and on the initial atomic state 
is examined. 

1 Introduction 
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& . The application of quantum optical devices in quantum information processing, 

specifically the efforts in constructing quantum memories [TJ [2J El S] and quan- 
tum networks [5] [6j [7] , is an active research area. The necessary part of such de- 
vices, employing the photons as flying q-bits, is the deterministic single-photon 
CN , emitter - the "photon pistol". The most promising scheme for the controlled 

generation of a single photon is based on the technique of vacuum stimulated 
Raman scattering involving adiabatic passage (STIRAP), proposed in [5] and 
then realized experimentally in [9] [TO] [IT] [12] . In these experiments the three- 
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t^J- , level A - type atom was adiabatically passing through a high-finesse cavity, one 

_-J ■ branch of the atomic A - type transitions was in resonance with the quantized 

f*^ , cavity field, while the other one was driven by the coherent laser field. In course 

CNj ■ of interaction with the fields inside the cavity the atom emitted a single photon. 

The two polarization degrees of freedom of the emitted photon provide the most 
natural way to encode the q-bit, however, in the early experiments [9] [10] the 
polarization of the emitted photons was not controlled, nor detected, and the 
rS ■ degeneracy of resonant levels was not taken into account. In the subsequent ex- 

C^ [ periments [TTJ [TJ] the magnetic field directed along the cavity axis was applied 

to lift the degeneracy of atomic levels and the frequencies of the driving field 
and of the cavity modes were adjusted to realize the non-degenerate A-schcmc 
with three Zeeman states, which allowed to produce the circularly polarized 
photons. The objective of this paper is to calculate the polarization state of the 
emitted photon and to study its dependence on the polarization of the driving 
laser field and on the atomic initial state for the arbitrary values of the level 
angular momenta and to look for the schemes enabling to produce photons with 
tailored polarizations. 

The STIRAP with the three-level atom with non-degenerate levels is well 
described in the reviews and textbooks (see, e.g., |T3l E]). In case of non- 
degenerate levels there exists the only dark state - the superposition of the 



ground state a and some metastable state b - uncoupled to the excited level c. 
In course of Raman scattering the atom is adiabatically transferred from the 
initial state a to the target state b. The STIRAP with degenerate atomic levels 
and with classical coherent resonant fields was studied in |15j , were it was shown 
that unlike the non-degenerate case the population from the initial ground state 
a is not always totally transferred to the target level 6, but only if the order 
of degeneracy of this target state is not less, than the order of degeneracy of 
the initial state. In our previous papers [TH1 E] we have studied the single- 
photon emission via non-adiabatic vacuum stimulated Raman scattering of short 
pulses on the atoms with degenerate levels and with totally or partially resolved 
hypcrfme structure, the polarization of the cavity mode was assumed to be 
well-determined and the dependence of the photon emission probability on the 
mutual orientation of polarizations of the driving pulse and that of the cavity 
mode was calculated. In the present paper the adiabatic Raman scattering 
(STIRAP) is considered, and it is assumed that the cavity equally sustains 
both polarization modes. 

In section 2 the interaction model and the evolution operator for this model 
in the adiabatic approximation are described, while in sections 3 and 4 the in- 
stantaneous eigenvectors of the interaction operator, which determine the evo- 
lution operator, are constructed and classified. In case of degenerate levels not 
only the number of these eigenvectors increases, but there appear the new types 
of these eigenvectors, non-existing in case of non-degenerate levels, like the dark 
states, which atomic part belongs to only one of the lower levels a or b, and 
the bright states, which couple the excited level c with only one of the lower 
levels a or b. In section 5 the formula for the polarization matrix of the emit- 
ted photon is obtained and in section 6 this formula is used for calculation of 
the photon polarization matrix for the transitions with the angular momenta 
J a = J c = 3, Jf, = 2 and Jb = J c = 1, J a = 2, corresponding to the transitions 
between the hyperfine structure components of the electronic levels 5Si/ 2 and 
5P 3 / 2 of the 85 Rb and s7 Rb atoms, which were employed in the experiments 
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2 Evolution operator 

Let us consider the coherent laser field with the carrier frequency u> c , which is 
in resonance with the frequency lo c q of an optically allowed transition J a — > J c 
between the ground state J a and the excited state J c , while the quantized cavity 
field with the carrier frequency u is in resonance with the frequency loq of an 
optically allowed transition Jb — s> J c between the long-lived state Jb and the 
same excited state J c (Fig.l). Here J a , Jb and J c are the values of the angular 
momenta of the levels. The coherent laser field is characterized by the electric 
field strength 

E c = e c (i)l c e-^*+c.c., (1) 
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Figure 1: The level diagram. 

while the quantized field of the cavity in the interaction representation is de- 
scribed by the operator: 



E = e(t)(oili + a 2 l 
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h.c. 



(2) 



where e c (t) and l c are the slowly varying amplitude and the unit polarization 
vector of the laser field, e(t) is the slowly varying amplitude of the cavity field, 
li and I2 are the two unit orthogonal vectors of the two polarization modes 
of this field, d\ and 0,2 are the photon annihilation operators for this modes. 
The temporal dependence e(t) of the cavity field amplitude appears due to 
the motion of the atom through the cavity or to some tailored alteration of 
cavity parameters. The equation for the slowly- varying density matrix p of the 
system, which consists of a single three-level atom and two-mode cavity field, 
in the rotating- wave approximation and in case of Raman resonance lu c o — lu c = 
luq — lu = A is as follows: 

d i 



dt P 2 



V(t),p 
V{t) = -2AP C + G{t) - 
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G(t) = n a (t)g a + n b {t)g b , 
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(3) 

(4) 
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Here P c is the projector on the subspace of the atomic excited level J c , Sl (t) = 
2\d a \e c (t)/Ti and ^(t) = 2\d b \e(t) / % are the reduced Rabi frequencies for the 
coherent laser field and for the cavity field, d a = d(J a J c ) and <4 = d(JbJ c ) 
are the reduced matrix elements of the electric dipole moment operator for the 
transitions J a — > J c and J b — » J c , while 
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g a and g & are the dimensionless electric dipole moment operators for the tran- 
sitions J a — > J c and Jb — > J c - These operators are expressed through Wigner 
3J-symbols and partial atomic operators 



Pm a m0 = \J a m a X Jpmp\, 


a, (3 


= a, b,c, 


in a following way [19]: 
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where l cq and l q are the circular components of polarization vectors l c and 1^. 
The solution of the equation ([3]) is expressed through the evolution operator 

S(t): 

p(t) = S(t)p(0)S + (t). (10) 

In the adiabatic approximation (see, e.g., [13]) the evolution operator S(t) is 
defined by the instantaneous eigenvectors \vk(t) > and eigenvalues Xk(t) of the 
interaction operator V(t): 

S{t) = Eexp{# fc (*)}M*) >< v k (0)\, (11) 



(t) = l f\ k {t')dt'. (12) 
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3 Bright states 

Since only the processes of a single photon emission are discussed in the present 
paper it is sufficient to limit the system space to the subspace with the basis 
vectors \J a m a > |0,0 >, \Jbmt, > |1,0 >, |J&m;, > 0, 1 >, |J c m c > |0,0 >, 
where \J a fn a >, | Jbfnb > and \J c m c > denote the atomic Zeeman states, while 
|ni,ri2 > (n-1,2 — 0, 1) are the field number states with ri\ photons in the first 
polarization mode and ri2 - in the second. This subspace with the dimension 
N = 2(J a + 2Jb + J c + 2) constitutes the invariant subspace of the interaction 
operator V(t), so that its matrix represents itself a square hermitian N x N 
matrix. The states \a > |0,0 > and \c > |0,0 >, which atomic part belongs to 
the level a or c, may be represented as columns with 2 J a + 1 or 2 J c + 1 elements 
correspondingly, while the states \b\ > |1, > +\b 2 > |0, 1 >, which atomic part 
belongs to the level b, may be represented as columns with 2(2 J b + 1) elements: 



|6i > |1,0>+|&2 > |0,1> = 



\bx> 

I&2 > 



Then the operator g a will be represented by the (2 J a + 1) x (2J C + 1) matrix, 
while the operator c/b will be represented by the 2(2Jf, + 1) x (2 J c + 1) matrix 

n - ( 9bi 

9b — 

\9b2 

were each block g^i (i = 1, 2) represents itself a (2 J& + 1) x (2 J c + 1) matrix. 

In order to find out the instantaneous eigenvectors \vk(t) > and eigenvalues 
Afc(t) of the interaction operator V(t) let us start with the eigenvectors of the 
operator G'(t)G(t), which acts at the subspace of the upper atomic level c. Let 
us denote as |D° > the states, if there are any, which are simultaneously the 
eigenvectors of both operators g\g a and g b gb with zero eigenvalues: 



.'/. 



{g a \D c n >=glg h \D c n >=Q, (13) 



the number of such sates being N%. These states remain uncoupled to the lower 
atomic levels a and b. Next, let us consider the states |C° >, which are the 
eigenvectors of operator g\g a with non-zero eigenvalues and at the same time 
the eigenvectors of the operator glgb with zero eigenvalues: 



fflffalC >= c\ n \C a n >, c an > 0, g\g b \C a n >= 0, (14) 

the number of such sates being iV° . These states are coupled to the lower atomic 
level a only. Similarly, the states |C„ >, coupled to the lower atomic level b 
only, are the eigenvectors of the operator glgb with non-zero eigenvalues, which 
are at the same time the eigenvectors of the operator glg a with zero eigenvalues: 

9 f b9b\C b n >= c 2 bn \C b n >, c bn > 0, glg a \C b n >= 0, (15) 

the number of such sates being N b . Finally, let us consider the states \C n (i) >, 
which satisfy the inequations: 

9a\C n {t) >^0, g b \C n (t)>^0, 

at any time. These states, coupled to both lower atomic levels a and b, may be 
obtained as the eigenvectors of the operator &(t)G(t) with non-zero eigenvalues 
at any time: 

G\t)G(t)\C n (t) >= c 2 n (t)\C n (t) >, c n (t) > 0, (16) 

the number of such sates being JV C . All the temporally independent N% states 
|.D£ >, N^ b states \C„ ,b > and temporally dependent N c states \C n (i) > con- 
stitute the complete (N% + N£ + N h c + N c = 2J C + 1) orthonormal set of states 
with the atomic part belonging to the atomic level c at any instant of time. 
The bright states 

\F% b >= -^-g a , b \C«> b >, n = 1, ..., N a c b , (17) 

C-a.bn 



\F n (t) >= J—G{t)\C n {t) >, n=l,...,N c , (18) 

coupled to the states \C^ b > and \C n (t) > by electric dipole transitions, also 
form an orthonormal set of states, as it follows from (jTT)) - ([T5]) . The atomic part 
of the temporally independent states \F*'° > belongs to the subspace of the 
only one lower level a or b, while the atomic part of the temporally dependent 
states \F n (t) > belongs to the subspace of both lower levels a and b. 

With the introduction of the states \D c n >, \C^ b >, \C n (t) > and \F%' b >, 
\F n (t) >, the interaction operator V(t) may be easily diagonalized. It has Nf? 
eigenvectors |D° >, then 2(iV ( ? + N b ) eigenvectors \V^ b ' n (t) >, which are linear 

superpositions of states \F^ b > and |C°' 6 >, and 2N C eigenvectors \Vn (t) > 
which are linear superpositions of states \F n (t) > and \C n (t) >. All these states 
\ D n >, \Vafil(t) > an d |K (±) (t) > constitute the orthonormal set of N f = 
2(2 J c + 1) — N* eigenvectors of the operator V(i) with non-zero eigenvalues. 



4 Dark states 

The other N d = N — N? eigenvectors \Dk(t) > of the interaction operator V(i) 
obtain zero eigenvalues: 

^(t)|U fc (t)>=0, k = l,...,N d . 

These states with the atomic part belonging to the subspace of the lower atomic 
levels a and b remain uncoupled to the upper level c. These states - dark states 
- satisfy the equation: 

G\t)\D k {t) >= 0, k = l,...,N d . (19) 

As it follows from (|T9j) . all the dark states \Dk(t) > are orthogonal to all the 
states \D c n >, |VJJ>(<) > and |V„ (±) (t) >. 

Among all the dark states let us distinguish first the dark states \D% >, 
which atomic part belongs to the lower level a only. These states are time 
independent and may be obtained as the eigenvectors of the operator g a g a with 
zero eigenvalues: 

g a gl\D%>=Q, fc=l,...,JV* (20) 

All the other eigenvectors \A n > of the operator g a g\ with non-zero eigenvalues, 
which satisfy the equation 

g a g\\A n >= a 2 n \A n >, a n > 0, n = 1, ..., ^V a , (21) 

constitute the subspace orthogonal to the subspace of the dark states |D£ >, 
its dimension being N a — 2J a + 1 — N%. Similarly the dark states \D b k >, 
which atomic part belongs to the lower level b only, may be obtained as the 
eigenvectors of the operator gtg b with zero eigenvalues: 

g b gl\D b k >=0,k = l,...,Nt (22) 



The N b = 2(2 Jb + 1) — N d eigenvectors \B n > of the operator gbg b with non-zero 
eigenvalues: 

g b gl\B n >= b 2 n \B n >, b n > 0, n = 1, ..., N b , (23) 

constitute the subspace orthogonal to the subspace of the dark states \D b k >. 

Let us now consider the dark states \D k b (t) >, which atomic part belongs 
to both atomic lower levels a and b. These states satisfy the equation 

&(t)\D?(t)>=0, k = l,...,N d b , (24) 

while 

§l\D?(t)>*0,fi\D?(t)>*0. 

The temporal dependence of these states may be immediately obtained from 
the equation (J24)) : 

\D?(t) >= z k (t)[n a (t)\B d k > -n b (t)\Ai >}, 

where Zk(t) is the normalization factor, while \Af > and \Bf, > are temporally 
independent states, which atomic parts belong to the levels a and b correspond- 
ingly and which satisfy the equation 

gt\A d k >=gi\B d k >^0. (25) 

Introducing the matrix 

N b 

D b = J2v\ B n><B n \, (26) 

n=l ™ 

containing only the eigenvectors \B n > of matrix gbg\ with non-zero eigenvalues, 
we may write the equation (J25I) as follows: 

\B d k >= D ba \Af >, b ba = D b g b gt (27) 

Now we may define the orthonormal set of states \A k > (k — 1, ..., N d b ) as the 
eigenvectors of the hermitian matrix Dl D ba with non-zero eigenvalues: 

Dl a Db a \A d k >=a 2 dk \A d k >,a dk >0. (28) 

Then, as it follows from ((27)) and (l28l) . the states a~^ k \B k > also constitute the 
orthonormal set, so that the orthonormal set of dark states \D k b (t) > may be 
expressed through the eigenvectors \A d > of the equation (|2"5|) by the formula: 

1 * ( } ^HT^PH ' (29) 



5 Single photon emission 

All the eigenvectors of the interaction operator V(t), comprising the set of states 

\ D n >» IK ( £(*) >> l^" (±) (0 > and the set of dark states \D a k ' b >, \Df (t) >, 
constitute the complete orthonormal set of states, which determines the evolu- 
tion operator (|11[) in the adiabatic approximation. Initially the atom is at its 
lower level a, while the cavity field is in its vacuum state, the initial atomic 
density matrix of the atom+field system being /5g . The classical coherent laser 
field is adiabatically switched on, while the interaction with the quantum field 
is adiabatically switched off in the time interval T, so that: 

o a (o) = o, n a {T) = n a0 , 
n b (o) = fibo, n 6 (T) = o. 

The atomic part of the initial states \Vn (0) > belongs to the atomic levels b 
and c only, so that only the states \Van (t) > and the dark states \D% > and 
\Df(t) > will contribute to the evolution operator ([11]). After the STIRAP 
process is finished at t = T the atomic part of the states | Van (T) > belongs 
to the levels a and c, the atomic part of the states \D% > belongs to the level 
a and the atomic part of the states \D^ b (T) > belongs to the level b, so that 
only the presence of the dark states |D£ 6 (i) > in the evolution operator con- 
tributes to the probability of the photon emission, while the presence of the 
states | Von {t) > and \D% > in the evolution operator reduces this probability. 
The total probability w of the single photon emission is 

w =< 1, 0\p F |1, > + < 0, 1|Af|0, 1 >, (30) 

while the polarization of the photon is described by the 2x2 polarization matrix 

l_ /<l,0|p F |l,0> <1,0|/3 F |0,1>\ . . 

a w V<0,l|y5 F |l,0> < 0, l\p F \0, 1 > J ' [ ' 

Here 

p F = tr A {Uffifrj, U = DbaP*, (32) 

the trace is carried out over the atomic variables, the operator D oa is defined 
by ([27]) and 

N d ab 

p d a =Y.—\ A i>< A iy ( 33 ) 

~. a-dk 
k=L 



where \Ai > and a^ are defined by 



6 Discussion 

Let the quantization axis Z be directed along the cavity axis and let us choose 
the two unit polarization vectors of the quantum field as the two unit vectors 



along the Cartesian axes X and Y: 

h = e x, I2 = e y, 

while the driving field propagates along the axis Y so that its polarization vector 
l c belongs to the plane XZ. In case of linearly polarized driving field 

l c = e z cos ip. + e x sin -0. 

The photon polarization matrix pip may be expressed through the Stokes pa- 
rameters £„ (n = 1, 2, 3) [20] : 

„ = 1 ( 1 + 6 & - *£» 
which determine the photon total degree of polarization 



In the experiments [9, 10 the resonant levels with angular momenta J a = J c = 3, 
Jb = 2 were the hyperfine structure components of the electronic levels 5S l 1 / 2 
and 5J3/2 of the S5 Rb atom, however the polarization of the driving field and 
the initial atomic state were not specified. So let the atom be initially in the 
equilibrium state 

fi = ^10,0X0,01, 

where P a is the projector on the subspace of the atomic ground state a. Then 
we obtain from the numeric calculations based on the formulae (|3"0"|) - ([3"3|) . that 
the photon emission probability w — 0.857 is less than unity and does not 
depend on the angle ip between the cavity axis Z and the polarization vector l c 
of the driving field, while the emitted photon is almost unpolarized, its degree 
of polarization is zero P = with the 7r-polarized driving field {ip = 0) and 
P = 0.07 with the a-polarized field (-0 = 7r/2). Now let us assume that the 
atom is initially prepared in a pure Zeeman state \J a = 3, m a = > with zero 
value of the angular momentum projection m a on the quantization axis Z, like 
in the experiment [TB]. Then we obtain from (J3"0")) - (|3"31 the dependence of the 
photon emission probability w on the angle tp presented at Figure 2. In this 
case with the 7r-polarized driving field (ip — 0) the photon is not emitted at all 
w = 0, while it is emitted with the unit probability w = 1 with the c-polarized 
field (if) = 7r/2) and at an angle ip — 0.685 rad (39°), at ip = 7r/2 its degree of 
polarization being P = 0.43 and at tp = 0.685 it is P = 0.80. 

In the experiments [HUH] the magnetic field directed along the cavity axis Z 
was applied to lift the degeneracy of atomic levels and the cavity field was tuned 
in resonance with the transition between the Zeeman state \F' = l,mp = > 
of the hyperfine structure component F' = 1 of the excited 5P3/2 electronic level 
of 87 Rb atom and one of the Zeeman states \F = 1, mp = ±1 > of the hyperfine 




angle \|/ (rad) 

Figure 2: The photon emission probability w on the transitions with the angular 
momenta J a = J c = 3, Jf, — 2 versus the angle "0 between the cavity axis Z and 
the polarization vector l c of the driving field in case of initial pure atomic state 
m a = with zero angular momentum projection on the quantization axis Z. 



structure component F — 1 of the ground level 551/2, while the linearly a- 
polarized driving held was in resonance with the transition between the same 
Zeeman state \F' — l,m,F> — > of the excited level and the other Zeeman 
state |.F = 1, mp = =pl > of the ground level, so that the emitted photons were 
circularly polarized, as it is evident from selection rules. 

In the experiment [18j the resonant levels with angular momenta Jf, = J c = 1, 
J a = 2 were the hyperhne structure components of the electronic levels 5Si/2 
and 5P3/2 of the 87 Rb atom and the atom was initially prepared in the pure 
Zeeman state \J a — 2, m a — >, while the driving held was linearly 7r-polarized. 
The results of numeric calculations of the photon emission probability w and the 
degree of polarization P versus the angle tp between the cavity axis Z and the 
polarization vector l c of the driving held for these transitions are presented at 
Figure 3. In case of 7r-polarized (-0 = 0) driving field, as it was in the experiment, 
the photon is emitted with unit probability w = 1, but it is fully unpolarized 
P = 0. The photon becomes unpolarized after the averaging over the states 
of the atom, however in the experiment |18) the same atom was driven twice 
in the time interval less than decoherence time to produce the sequence of two 
photons with strongly correlated polarizations. As it may be seen from Figure 3, 
with the same experimental setup but with the driving field linearly er-polarized 
(tp = 7r/2) the emitted photon will be fully polarized P = 1 (linearly polarized 
along the axis Y as follows from calculations), though the emission probability 
in this case is reduced to w = 0.25. To obtain the linearly polarized photons 
with unit probability the scheme with lower values of angular momenta J a = 0, 
Jb = Jc = 1, will be suitable. In this case the cr-polarized driving field will 
produce the photons linearly polarized along the axis Y with unit probability. 
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Figure 3: The photon emission probability w (solid line) and the degree of 
polarization P (dashed line) on the transitions with the angular momenta Jb = 
J c = 1, J a = 2 versus the angle ip between the cavity axis Z and the polarization 
vector l c of the driving field in case of initial pure atomic state m a = with 
zero angular momentum projection on the quantization axis Z. 

7 Conclusions 

In the present paper the general formulae for calculation of the polarization 
matrix of a single photon emitted in the microcavity by a single three-level A- 
type atom with degenerate levels, driven adiabatically by the classical laser field 
from the ground state a to the long-lived state b through the excited state c, are 
obtained for the arbitrary polarization of the driving field and arbitrary values 
of the angular momenta J a , Jb and J c of the resonant atomic levels. 

For the transitions with the angular momenta J a = J c = 3, Jf, = 2 and 
Jb = J c = 1, J a = 2, corresponding to the transitions between the hyperfine 
structure components of the electronic levels 5Si / 2 and 5P 3 / 2 of the S5 Rb and 
87 Rb atoms, which were employed in the experiments [HI [TOJ [TS], and for the 
atom initially prepared at pure Zeeman state with zero projection m a = on 
the quantization axis, the dependencies of the photon emission probability w 
and degree of polarization P on the angle ip between the cavity axis and the 
polarization vector l c of the driving field were calculated numerically. In case 
of transitions with J a = J c = 3, Jb = 2, the unit emission probability may be 
obtained at the angles ip = 39° and ip = 90° with the degrees of polarization 
P = 0.80 and P — 0.43 correspondingly. In case of transitions with Jb = J c = 1, 
J a = 2, at tp = the fully unpolarized photon is emitted with unit probability, 
while at tp = 90° fully linearly polarized photon is emitted with the probability 
w = 0.25. 
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The fully linearly polarized photon emitted with unit probability may be 
obtained by means of cr-polarized driving field (ip = 90°) with lower values of 
the level angular momenta J(, = J c = 1, J a = 0. 
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